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Abstract 

Bifurcation cascades in conservative systems are shown to exhibit a generalized diagram, which contains all relevant 
informations regarding the location of periodic orbits (resonances), their width (island size), irrational tori and the 
infinite higher-order resonances, showing the intricate way they are born. Contraction rates for islands sizes, along 
period-doubling bifurcations, are estimated to be ai ~ 3.9. Results are demonstrated for the standard map and for the 
continuous Henon-Heiles potential. The methods used here are very suitable to find periodic orbits in conservative 
systems, and to characterize the regular, mixed or chaotic dynamics as the nonlinear parameter is varied. 
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1. Introduction 

Almost all physical systems in nature are so com- 
plex that long time predictions are nigh impossible. 
This is a characteristic of nonintegrable chaotic systems 
whose eff'ects are visible in celestial mechanics, plasma 
physics, general relativity, quantum physics, commu- 
nications problems, hard beats, social and stock mar- 
ket behaviors, weather forecast, among others. Thus, 
the precise description of the dynamics in nonintegrable 
systems is essential for the understanding of nature. 
One intrinsic and fundamental phenomenum of such 
nonintegrable systems, is that realistic stable orbits may 
vanish (or be born) when the nonlinear parameter varies. 
This can lead to a complicated behavior with a cas- 
cade of new orbits, which is satisfactory described by 
a bifurcation diagram [1]. One example is the period 
doubling bifurcation (PDB) cascade, well understood 
in one dimensional dissipative discrete systems con- 
taining one parameter, where the intervals in the pa- 
rameter, between successive PDBs, tend to a geomet- 
ric progression with an universal ratio of 1 = 1 /4.66, 
with d being the Feigenbaum constant. For two dimen- 
sional dissipative systems, bifurcation cascades mani- 
fest themselves inside periodic stable structures in the 
two parameter space, which appear to be generic for a 
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large class of systems. Shrimp-like structures are one 
example, and appear in the dissipatives Henon [2, 3] 
and ratchet [4] discrete systems, in continuous mod- 
els [5], among many others. Such structures allow an 
analysis of geometric approximation ratios, and pro- 
vide a very clear understanding of the dissipative dy- 
namics. In contrast to dissipative systems, for conser- 
vative nonintegrable systems the description of bifurca- 
tions cascades is much more complicated. It is known 
[1, 6], that by magnifying stable points in the phase 
space, a mixture of surrounding stable and unstable 
fixed points is found. This repeats itself for every stable 
fixed point, as explained by the Poincare-Birkhoff* theo- 
rem [7] . How higher order stable orbits bifurcate, the 
islands around them vary and are interconnected, and 
irrational tori behave as the nonlinear parameter of the 
system changes, is an interesting problem that still de- 
serves to be deeply investigated. For 2 degrees of free- 
dom area-preserving discrete systems, the intervals in 
the parameter between successive PDBs, tend to a ge- 
ometric progression with a ratio of \I6h ~ 1/8.72 [1]. 
Results have also been extended to higher-dimensional 
systems (please see [8, 9, 10, 11] for more details). 

This work uses convergence properties of the Finite 
Time Lyapunov Exponent (FTLE) to explore the dy- 
namics of conservative systems in a mixed plot: initial 
condition (IC) versus the nonlinear parameter. The lo- 
cation of stable orbital points is easily found in such 
plot, and conservative bifurcation diagrams recognized 
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to have a generalized form, containing infinite sub- 
diagrams with all rational/irrational tori from the peri- 
odic orbits (POs), independent of the period. This is a 
nice complete description, and extension, from an early 
work [12], done for another dynamical system. Re- 
sults are remarkable and show the very complex, self- 
similar and generic bifurcation structure in conserva- 
tive systems with only one parameter. They also sug- 
gest that contraction rates for the islands sizes, along 
PDBs, approach the constant aj ~ 3.9. A detailed nu- 
merical analysis is performed for the standard map, in- 
cluding the use of the first recurrence times instead the 
FTLEs, to demonstrate that results are independent of 
the method. A generalization is shown for the continu- 
ous Henon-Heiles potential. 

2. Discrete Model: The Standard Map 

To start it is appropriate to present results using a 
well known general model with wide applications, the 
Chirikov standard mapping, which is given by [13]: 

f Pn+i= Pn-^ (K/271) sin(27rx„) mod 1, 

\ Xn+l = Xn+ Pn+l mod 1. 

K is the nonlinearity parameter and are respec- 

tively position and momentum at discrete times n = 
1,2, ...,N. It is known that period- 1 (shortly writ- 
ten per-1) fixed points slyq pi = l/2m (m integer) and 
xi = 0, ±1/2. The point xi = is always unstable while 
xi = ±1/2 becomes unstable for K > 4. There exist 
also per-1 fixed points related to accelerator modes [13] 
whose stability condition is |2 ± ^cosxi/| < 2, with 
K sin xii = 2nl and / integer. For higher periods there 
are primary families of periodic points (which exist in 
the limit ^ ^ 0) and bifurcation families which are 
born only for larger values of K (See [1] for more de- 
tails). 



torus close to xq. Assuming that (Axq)^ ^ 0.0, the stan- 
dard map can be linearized around x^ and the FTLE de- 
termined analytically, after n iterations, from the eigen- 
values of the composed Jacobian {(1, -K), (1,1- K)}^. 
This FTLE is plotted in Fig. 1(a) as a black continuous 
line, and it converges exactly to zero only when n ^ oo. 
This means that ICs from stable tori around the fixed 
point, take a longer time to converge exactly to zero 
than the fixed point itself. This behavior can also be 
observed numerically by determining the FTLEs using 
Benettin's algorithm [14, 15], which includes the Gram- 
Schmidt re-orthonormalization procedure. We use the 
six exemplary orbits shown in the inset of Fig. 1(a): the 
stable fixed point, demarked with a cross, and the five 
irrational tori around the fixed point. All trajectories 
are regular and have LEs exactly equal zero for infinite 
times. However, when calculating the FTLEs for the 
distinct ICs, we observe that ICs closer to the fixed point 
converge faster to zero. In Fig. 1(a), the decay curves for 
the FTLEs versus times are plotted (starting from be- 
low) for the six ICs shown in the inset (starting from the 
fixed point). As ICs are taken more away from the pe- 
riodic point, FTLEs converge slower to zero. The mag- 
nitude of the FTLEs between diff'erent irrational tori are 
really small and not significant for any purpose. How- 
ever, for the IC exactly on the PO, the FTLE converges 
faster to zero than other ICs around it. Essential to men- 
tion is that this is not a numerical convergence artefact 
due to the numerical method, but an analytical property, 
as shown above (continuous line in Fig. 1(a)). Thus, 




3. Method: Finite Time Lyapunov Exponents 

The key idea for the success of our proposal is the ob- 
servation that the numerical convergence of the FTLEs 
is distinct for diff'erent ICs (same nonlinear parameter 
K), even between regular trajectories. To make this 
clear, consider, for example, that the IC is exactly on the 
stable (for ^ < 4) fixed point xq = 1/2. It can be calcu- 
lated analytically that the corresponding Lyapunov ex- 
ponent (LE), after one iteration, is exactly zero. But now 
consider that we start with an IC close to this fixed point, 
say Xq = 1/2 + Axq, which can be an irrational regular 



Figure 1 : (Color online) (a) Log-log plot of the FTLEs (starting from 
below) as a function of the iteration time for the six trajectories shown 
in the inset (starting from the fixed point), and (b) phase space for 
K = 2.6. 

even though such small FTLEs are insignificant to dis- 
tinguish between the irrational tori, they allow us to rec- 
ognize where POs are located in phase space, and to 
understand the very complex and self- similar behaviors, 
which occur close to the POs as the nonlinear param- 
eter K changes. For larger times, FTLEs in Fig. 1(a) 
continue to decrease linearly to zero, until the machine 
precision is reached, and when the actual method can- 
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not be used anymore. Even though the method was 
explained using the standard map, it is equally applied 
to other conservative dynamical systems. 

4. The Generalized Diagrams 

Using the properties explained in Section 3, a very 
clarifying plot can be constructed, which allows us to 
recognize the bifurcation diagram in conservative sys- 
tems in a simple way. Figure 2 shows the FTLE (see 
colors bar) in the mixed space xq x K with po = 0.0. 
Black lines are related to those ICs for which the FTLEs 
converge faster to zero and POs exist (this was checked 
explicitly for many black lines). Dark to light yellow 
points, around the main black lines, are related to irra- 
tional tori and also define the size of the correspond- 
ing island. These are regular trajectories for which the 
associated FTLEs are still not that close to zero, as 
those FTLEs from ICs which start exactly on the black 
lines (POs). For n ^ oo these points also become 
black. Light yellow, blue to red points, are ICs related 
to chaotic trajectories. Clearly we observe that the fixed 
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Figure 2: (Color online) FTLEs in the mixed space xq x/T for the stan- 
dard map with po = 0.0, a grid of 10^ x 10^ points and 10^ iterations. 

point xq = xi = 0.0 is unstable (positive FTLEs along 
the line xq = 0.0) for any K value. In the chaotic re- 
gion, FTLEs increase monotonically with K and the dy- 
namics is not much exciting. However, in the regular 
to chaotic transition region, trajectories combine them- 
selves in a very complicated way, bifurcations occur and 
the dynamics is very rich as K changes. For clarifica- 
tion. Fig. 1(b) shows the corresponding phase space dy- 
namics for K = 2.6 [see white dashed line in Fig. 2]. 
One per-1 (at x = ±0.5, /?o = 0-0) and one per-4 (see 
numbers i, ii, iii, iv showing the orbital points) POs are 
shown. The per-1 orbital points correspond to the black 
points in Fig. 2 located at xq = ±0.5 and K = 2.6. This 
is the primary resonance. The two orbital points (ii,iii), 
from the per-4 PO, belong to the black curves in Fig. 2. 



Since we always started along the initial line po = 0.0, 
Fig. 2 is able to detect only the two orbital points (ii,iii) 
along this line. The whole dark to yellow region around 
xo ± 0.5 (for ^ < 4) in Fig. 2, gives the island size and 
is related to all irrational tori around this PO. 

ICs which leave to stable POs (black lines), are char- 
acterized by the "Christmas tree" structure around it. 
For example, for K = 2.0 in Fig. 2, two per-4 trees 
are born, one at xq = -0.5, and the other one at 
xq = 0.5. All trees have a main "stem" in the center, 
which is black and gives the location of the stable PO 
when K varies. Since both mentioned trees are born 
from Xq = ±0.5, they are secondary resonances with 
per-4. A sequence of identical trees (secondary reso- 
nances) with smaller and smaller sizes, and increasing 
periods 8, 12, 16, . . . [see numbers on the sequence of 
trees Fig. 2], are observed by decreasing K. The stem 
of these trees, which are all born from the PO x = ±0.5, 
are related to the rational tori. The width of the tree is 
the width of the island along the line po = 0.0. In addi- 
tion, these trees approach to each other as they get closer 
to the unstable fixed point at xq = 0.0 and K 0.0. A 
large sequence of approximating trees is also observed 
(not shown here) close to ^ ~ 4.0 (where the per-1 PO 
becomes unstable) and xq ±0.5. The main stems 
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Figure 3: Magnifications of (a) the per-4 fork from the box from Fig. 2 
and (b) per- 8 fork from box from (a). Numbers indicate the periods of 
the POs. 

in the above trees exist until they become unstable at 
the very end of the tree (on the right), where a pitch- 
fork bifurcation occurs. This is better visualized in the 
magnification shown in Fig. 3(a) which, roughly speak- 
ing, summarizes results contained in Figs. 7 and 8 from 
[12], obtained for another dynamical system. However, 
Fig. 3(a) contains many essential additional informa- 
tions about the rich underline dynamics, as described 
below. The main stem (stable PO) becomes unstable 
close to ^ ~ 3.145 and xq ~ 0.25, as can be recog- 
nized by the yellow points (see dashed line). Simultane- 
ously around this point, two black lines (stable POs) are 
born with the same per-4, characterizing the pitchfork 
bifurcation. In fact, a "fork"-like structure (diagram) is 
born at the bifurcation point and contains two new main 
trees. This fork is the generalized diagram, which lives 
in the entire portion of the mixed plot where the tran- 
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sition from regular to chaotic motion occurs. This was 
checked for many other magnifications, even for higher 
order periods. Since the interesting dynamics and all bi- 
furcation processes occur through the fork and it trees, 
it is enough to discuss just one fork to understand the 
very rich dynamics from the whole Fig. 2. 

For this we start discussing the magnification of the 
per-4 fork shown in Fig. 3(a). Inside one fork, two 
main trees always exist: one called the Line Tree (LT), 
where PDBs occur for ICs on the line po = 0.0; and 
the other is the Outline Tree (OT), where PDBs oc- 
cur for ICs outside po = 0.0. In the example from 
Fig. 3(a), the OT is the lower one and the LT is the up- 
per one, which is identical to all secondary trees shown 
in Fig. 2. LT and OT have a main stem related to dis- 
tinct POs, but with the same period [see the numbers 
in Fig. 3(a)]. Emerging almost perpendicularly from 
the main stems, infinite smaller sub-trees are born be- 
low and above. Such new sub-trees have also a main 
stem related to another PO with higher period. In the 
ramifications shown in Fig. 3(a), the new POs emerg- 
ing from the LT, as K decreases, are: 8 (new fork, see 
box), 16, 24, 32, . . . while for the OT are: 20, 28, 36, . . .. 
These new ramifications (sub-trees) are the rational tori 
(higher-order resonances) which live around the per-4 
point. At the end of the fork (on the right), the main 
stems from the LT and OT suff'er a PDB phenomenum. 
This is shown in details in Fig. 3(b), which is a mag- 
nification of the white box per- 8 fork from Fig. 3(a). 
Besides the new trees emerging from the main stem, a 
continuous transition to lighter colors occurs when we 
go away from the main stem, in the direction of the 
border of the tree. All these points are related to ir- 
rational tori around the main PO. At the border of all 
trees, the breakdown of the last irrational torus, of the 
corresponding PO, can be observed by the transition 
from yellow to red points, where positive FTLEs ap- 
pear and the dynamics becomes chaotic. Remarkable is 
that the dynamics around this new unstable point, as K 
increases, is similar to the dynamics around the unsta- 
ble per-1 (xq = 0.0) orbit, observed in Fig. 2 for K close 
to zero. Collecting the main properties of the general- 
ized fork, suppose it has a per-^ for the main stems of 
the LT and OT. The ramifications have distinct periods, 
namely for the LT we have, starting from the right bor- 
der of the fork: 2q (new fork), and the ramifications re- 
lated to the rational tori have periods 4q, 6q, 8^, . . .. For 
the OT we have 2q (not seen along the line po = 0.0) 
with the ramifications 5q, 7q, 9q, . . .. From successive 
magnifications (not shown) of the generalized fork, we 
found that it contracts with a ratio Sp ~ 8.7, and that 
the islands sizes, at the PDBs, approach the estimated 



contraction rate aj = dq/d2q ~ 3.90 ± 0.04, where dq 
is the island size for ^ = 128. In order to estimate 
dq, we magnified the forks with main periods ^ = 128 
and 256, using times n = 10' . For such times, the bor- 
der line between the last irrational tori around the POs 
and the chaotic motion becomes evident, and dq can be 
determined from the numerical data of the FTLEs. 

5. A continuous model 

The general statement of our results can be confirmed 
showing the existence of trees in the continuous Henon- 
Heiles (HH) potential V(x,y) = (x^-\-y^)/2-\-x^y-y^ /3. It 
is a truncation of the Toda potential [1] used in celestial 
mechanics and extensively mentioned as a benchmark 
to nonlinear conservative dynamics. In this case the 
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Figure 4: FTLE in the space E xyo for the HH potential with xq = 
= 0.0, while is obtained from energy conservation. 

nonlinear parameter is the total energy E. For E = 1/24 
the dynamics is regular, and when E ^ 1 / 6 it becomes 
mixed and finally chaotic. Figure 4 shows the FTLEs 
in a portion of the the mixed plot E x yo for the HH 
problem. As in Fig. 2, black lines are related to those 
ICs which leave to close to zero FTLEs, where POs are 
present. Dark to yellow colors are related to the ratio- 
nal/irrational tori. Clearly we observe the existence of 
the stable trees, the forks (see white box for one ex- 
ample) and the infinite sequence of trees ramifications 
(higher order resonances), always appearing in form of 
the tree described in details for the standard map. 

6. Changing the method: First recurrence times 

The purpose of this section is to show that the gener- 
alized diagrams can be obtained using another method. 
Despite results using the FTLEs are totally trustworthy, 
it is proper to use a distinct method. For this we repro- 
duce Fig. 2 using the first recurrence times (FRTs) in- 
stead the FTLEs. We start the simulation with ICs inside 
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the interval -0.5 < xq < 0.5, po = 0.0, iterate the stan- 
dard map (1) and reckon Axn = xq - Xn, Apn = Po - Pn- 
When |Ax„| < 10""^ and \Apn\ < 10""^, we record Ti = n. 
For n > Tfnax = 10"^ we stop the simulation. In this 
way we compute the number of iterations Ti , which a 
given trajectory needs to come back to the initial point 
(inside a circle of radius 10""^) for the first time. It is 
expected that for POs the time Ti is smaller (equal to 
the period of the orbit), and that for chaotic trajectories 
the limit n = T^ax is reached. In between, all values of 
Ti are possible, and are usually related to irrational tori 
or POs with very large periods. Results are shown in 
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Figure 5: FRT Ti (see color bar) for the standard map (1) in the plot 
ICs versus K. 

Fig. 5. Plotted is Ti (see color bar) as a function of ICs 
xo and K. Black points are related to smaller times and 
are identified as the POs. Purple to red points are the 
irrational tori with intermediate values of Ti. Yellow 
points are related to chaotic trajectories, which never 
returned (for n < 10^) to the IC. Even though Fig. 5 is 
plotted with diff'erent colors, all trees, forks structures, 
widths and scalings, are identical to Fig. 2. In addition, 
the FRTs are able to find the location of unstable POs, 
as for example, the straight horizontal line at x = 0.0, 
which is the unstable fixed point for any K. 

The reason why the FRTs work so well to reproduce 
the results from the FTLEs is a discussion which we 
cannot explore in a satisfactory way in this work. It has 
also other important properties. The main point here 
is to demonstrate that our results are real physical ef- 
fects, relevant and independent on the method: FTLEs, 
FRTs, or any other quantity whose time convergence is 
aff'ected by the ICs. 

7. Conclusions 

To conclude, period doubling bifurcations in conser- 
vative systems are shown to exhibit a generalized stable 
diagram in the plot: initial conditions versus the nonlin- 
ear parameter. This diagram is composed of two main 



trees whose main stems are the periodic orbits (reso- 
nances) and infinite ramifications, which are all ratio- 
nal tori (higher-order resonances) surrounding the pe- 
riodic orbits. Results are discussed in details for the 
standard map and extended to the continuous Henon- 
Heiles potential. The generalized trees were also con- 
firmed to exist (not shown), for the nontwist standard 
map [16] and the hydrogen atom subjected to an exter- 
nal magnetic field [17], and can be recognized in [12], 
giving additional evidences about their wide applicabil- 
ity in conservative systems. First recurrence times and 
finite time Lyapunov exponents were used for the stan- 
dard map, showing that results are independent of the 
method. Other methods could be used to find the trees 
as, for example, the finite time rotation numbers. This 
is what Fig. 5, from the recent work about the nontwist 
standard map [18], suggests, providing one more exam- 
ple of a distinct system, with another method, where 
the generalized diagram should appear. Our results are 
also independent of the periods of the periodic orbits, 
and provide general clues for conservative systems con- 
cerning: the regular/mixed/chaotic dynamics, the loca- 
tion of stable orbital points, the generic mechanism of 
bifurcation phenomenum, and that contraction rates of 
islands sizes, along period doubling bifurcations, fol- 
low a I ~ 3.9. To finish, from the present analysis it 
is possible to conjecture that, for higher dimensional 
systems, the dynamics outside the generalized bifurca- 
tion diagram, and for small nonlinearities, should pro- 
duce Arnold stripes [19], which display the ICs from 
the stochastic Arnold web as a function of the nolinear- 
ity parameter. 
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